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ABSTRACT 


A mathematical model has been developed for simulation 
of forced circulation evaporation process. Computer programme 
written in FORTRAN, based on the mathematical model predicts 
the performance of forced circulation evaporator for variation . 
in operational parameters like feed rate, feed temperature, 
concentration of feed liquor and recycle ratio. Mathematical 
model jjredicts various quantities like temperature, pressure, 
pressure gradient, tube side film coefficient, overall heat 
transfer coefficient, liquor concentration, void fraction at 
each point along the tube length. Overall performance is 
determined by integrating the longitudinal profiles. Mathe- 
matical model developed is based on basic continuity, momentum, 
and energy balance equations applied to differential segments 
of the tube. Along with basic balance equations, different 
correlations for pressure drop and heat transfer coefficients 
for the non-boiling and boiling zones have been used in the 
model. Correlations for engineering properties of black 
liquor, considered as process liquor, are Available from earlier 
work. 

Due to lack of suitable plant data, model was tested 
on Badger's experimental data on vaporization of water in 
long vertical tube. The model was found to predict Badger's 
result satisfactorily. Model was used to study the performance 
of forced circulation evaporation in a single vertical tube, 
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and in a three pass unit with single tube per pass for both 
vertical and horizontal orientation of calendria. Bamboo 
and Bagasse black liquors were considered to illustrate 
the effect of viscosity on the performance of forced circu- 
lation evaporator. 
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CHAPTER 1 


INTRODUCTION 

Evaporation is an important chemical engineering 
operation used for concentrating solutions in diverse appli- 
cations in chemical, food, pharmaceutical and related industries. 
A summary of the general features of some of the more popular 
types of evaporators is given in Table 1.1. Forced circulation 
evaporators (FCE) find applications in concentrating moderately 
viscous solutions (80-100 cp) with tendencies for solid 
separation and/or scale formation. Typical process applications 
of FCE include evaporation of salt brines, electrolytic caustic 
solutions, black liquor, aluminate liquor etc. 

The heating element (calendria) of a FCE is essentially 
a multipass (2-6 passes) shell and tube heat exchanger and 
the heated process liquor is flashed in a separator. The 
desired increase in concentration of the liquors is obtained 
by recirculating the liquor after flashing. Liquor velocity in 
the calendria tubes is maintained at 1-6 m/s with a centri- 
fugal pump. Calendria can be positioned vertically or hori- 
zontally depending upon space availability and maintenance 
considerations. Figure 1.1 and Figure 1.2 [27 ] give schematic 
sketches of FCE with vertical and horizontal calendria, 
respectively. 
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TABLE 1.1 : GENERAL FEATURES OE EVAPORATORS USED IH 
CHEMICAL PROCESS INDUSTRIES 


Peature 

SE 

LTVE 

PPE 

r FCE 

ATPE 

Orientation 

Yert. 

Yert. 

Yert 

Vert/Hor. 

Vert/Hor. 

Tube length, 
m 

1-2 

6-12 

2-10 

2-7 

2-6 

Tube dia. ,mm 

25-50 

25-50 

25-50 

25-50 

300-2000 

Liquor Recir- 
culation 

HO 

HO 

HO 

YES 

HO 

Limiting 
viscosity ,cp 

50 

15 

50 

80-100 

3000 p 

U, value 
kW/m 2 °K 0 

.8-2.8 

1-3.4 

0. 6-1.7 

2-11 

Variable 8, 

Applications: 
scaling liq. 

■ HO 

YES 

HO 

YES 

YES 

Heat sensitive 
liquor 

HO 

HO 

YES 

HO 

YES 

Foaming liquor 

HO 

YES 

HO 

YES 

YES 

Slurrie s 

HO 

HO 

HO 

YES 

YES 


SE - standard evaporator 
LTVE - long tube vertical evaporator 
EPE - falling film evaporator 
PCE - forced circulation evaporator 
ATPE - agitated thin film evaporator 

a Strongly dependent upon viscosity 




Dual density 
mesh mist 
extractor 


Impingement 

baffle 


Steam 


Vortex I 
breaker 


Condensate 


Liquid concentrate from this line (by inserting a 


Forced-circulation vertical 
evaporator. * 
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Simulation of an evaporator unit consists in predicting 
profiles of liquor temperature, concentration, heat flux, 
overall heat transfer coefficient and evaporation rate* In 
simulating the performance of an FOB unit process conditions 
such as flowrate, temperature, concentration, inlet pressure 
of the feed liquor and vacuum in the seperator, recycle ratio, 
and steam pressure are generally available. A computer 
programme is developed in this work to predict the effect of 
variations in concentration, temperature of feed liquor and 
recycle ratio on the performance of a BCE unit. Mathematical 
model required for the simulation is obtained for steady 
state operation of the evaporator. The model is developed 
using basic continuity, momentum, and energy balance equations, 
correlations for pressure drop, hold-up and heat transfer in 
the non-boiling and boiling zones of the evaporator and the 
engineering properties of the process liquor. The program 
gave satisfactory predictions of Badger's [10 ] boiling heat 
transfer results with minor modifications to suit experimental 
conditions. The model was applied to a single tube vertical 
BCE to study the effect on performance of feed temperature 
(90-110°C), concentration (35-60 per cent) and rate (5-10 kg/s) 
for the concentration of black liquor. The model was also 
used to study the effects of feed concentrations, feed 
temperature, and recycle ratio, on the performance of a three 
pass FCE for both horizontal and vertical calendria arrangements 
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for processing black liquor. The effect of viscosity on 
performance of FCE was illustrated by bagasse black liquor 
which has a higher viscosity compared to bamboo black 
liquor. 



CHAPTER 2 


HYDRODYNAMICS AIR) HEAT TRANSFER IN FORCED CIRCULATION EVAPORATOR 

Performance of an evaporator depends on temperature 
drop and overall heat transfer coefficient. Temperature 
drop is determined by the temperatures of steam and the 
boiling liquid and boiling point rise. Overall heat transfer 
coefficient is strongly influenced by equipment geometry, 
hydrodynamics of fluid flow and various liquid properties 
like viscosity, thermal conductivity, specific heat, and 
operation of the unit. During evaporation of viscous solutions 
over the high concentration range, the resistance of the 
liquor side film becomes a major factor controlling the overall 
he at.. transfer coefficient. 

2.1 Porced Circulation (PC) and Natural Circulation (NC) Units ; 

NC and PC units form two complementary systems of 
evaporators in many process applications. The ' former generally 
constitutes the individual effects of a multiple effect 
evaporation plant, providing a large increase in liquor 
concentration. Over the high concentration range the liquors 
are generally viscous (above 10 cp) and give low heat transfer 
coefficients in NC units. This limitation is overcome by a 
PC unit, which circulates the liquor through the tubes at a 
higher velocity (1-6 m/s) and improves the rate of heat 
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transfer. _ For example, in a pulp mill, the evaporation 
plant consists of * a system of multiple effect evaporators having 
4-6 effects for the concentration range of 15-55 per cent 
solids and a FC finisher unit for the range of 55-65 per cent 
solids. Agarwal [ 9 ] has developed a simulation program for 
a 1TC long tube vertical evaporator (LTV) typical of the former. 
This investigation deals primarily with a model for FC eva- 
porator. 

Liquor flow pattern in a NC, LTV unit consists of a 
homogeneous solution in the non-boiling zone and two-phase 
vapor liquid mixture in the boiling .. zone . Flow pattern in the 
boiling zone consists of bubble flow, slug flow and annular 
flow regions. Agarwal has developed a model based on basic 
momentum, mass and energy balance equations and correlations 
describing pressure drop and heat transfer characteristics. 
Lochhart-Martinelli [16] correlation was used for determining 
frictional pressure drop for the two-phase mixture in the 
boiling zone. Coulson-McNelly ' s [ll] and Penman-Tait 1 s [14] 
correlations were used for the nucleate boiling and annular 
boiling regions respectively. A maximum in the longitudinal 
temperature profile was used as the criteria for the transition 
from non-boiling to boiling zone. A minimum in the profile of 
pressure gradient [ 13 ] was used to determine the transition 
between nucleate boiling and annular film boiling regions. 

The model adequately represented the boiling heat transfer 
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data of Badger [10 ] and also experimental longitudinal 
profiles of temperature and heat transfer coefficient of 
Gudmundson [19 ] . fluid flow and heat transfer mechanisms of 
FCE units depend upon the vertical or horizontal orientation 
of the calendria and are also influenced by the liquor 
viscosity, and velocity. These mechanisms are discussed in 
the next section. 

2.2 Hydrodynamics and Heat Transfer ; 

2,2.1 Plow Pattern in FCE: 

Heat transfer in evaporators is intimately related 
to fluid flow pattern. In the non-boiling zone of FCE fluid 
is a homogeneous liquid generally in turbulent flow except 
while handling viscous solutions when the flow can become viscous. 
In the boiling zone fluid flow pattern is different from that 
in non-boiling zone. In the boiling zone simultaneous flow 
of vapour and liquid phases occurs, and the inherent pressure 
change along the tube continuously changes the phase distri- 
bution to give varying longitudinal flow patterns, with 
corresponding changes in heat transfer mechanisms. A complete 
hydrodynamic description of the two-phase flow pattern 
requires a knowledge of the variation of void fraction and 
distribution of velocity and shear forces in the vapor and 
liquid phases of the system. Flow patterns in a two -phase 
vapor-liquid system have been classified according to visual 
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observation [1.5,1] and also empirically by Lockhart, and Martinell 
[l6 ] by a combination of laminar or turbulent flow conditions, 
prevailing for the superficial flow of either phase in the 
conduit. The spectrum of visually observed flow patterns for 
two-phase vapor liquid flow is described below. 

Bubble Flow : Bubble flow pattern occurs at low void fraction. 

The liquid phase is continuous and vapor phase distributed 
in the form of bubbles forming the discontinuous phase in 
this regime, followed by nucleate boiling heat transfer 
mechanism. 

S lug Flow s In slug flow ‘large slugs of vapor appear in the 
flow as a result of agglomeration of vapor bubbles. The 
slugs of vapor accelerate faster than compared to liquid flow 
and create slip between phases. Slug of liquid and vapor 
alternate each other and gives rise to instability. This 
flow pattern would occur at moderate void fractions and 
relatively low flow velocities. Heat transfer occurs both by 
nucleate boiling and forced convective mechanisms. 

Annular Flow ; Increase in vapor velocity results in a central 

core of vapor flow and an annular climbing film of liquor 

and heat transfer across the liquid occurs by forced convection. 

Mist Blow ; With an increase in temperature gradient dry 
patches begin to appear on the tube surface and liquor 
droplets are carried as mist by high velocity vapor. 
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In the boiling zone heat transfer mechanisms are not 
very distinct, and there is considerable overlap between 
various mechanisms in the boiling zone. Heat transfer 
occurs both by nucleate boiling mechanism and forced 
convection. The convective heat transfer in the boiling 
zone does not prevent the wall temperature from rising 
above saturation temperature [21]. As vaporization 
increases additional nucleation sites are activated and the 
number of nucleation sites increase and fully developed 
nucleate boiling occurs. In FC evaporators with vertical 
calendria the flow pattern is limited to bubble flow as vapor 
bubbles have little chance to coalesce to form slugs. 

Coexistence of nucleate boiling and forced convection heat 
transfer in forced convection vaporization has been experi- 
mentally verified [28]* Flow pattern in FCE with horizontal 
calendria consists of turbulent flow in non-boiling zone, 
and bubble flow in boiling zone with forced convection and 
nucleate boiling as the dominant heat transfer mechanisms in 
the boiling zone. 

lockhart-Martinelli classified the two-phase flow 
patterns in boiling zone empirically by a combination of 
laminar or turbulent flow conditions prevailing for the 
superficial flow of either phase in the conduit. Consequently 
these flow mechanisms result in laminar -laminar , laminar- 
turbulent, turbulent-laminar and turbulent -turbulent conditions. 
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Lockhart-Martinelli parameter is used in beet transfer 
correlations in Section 2.2.3 for the boiling section. 


2.2.2 Pressure Drop Correlations ; 

Pressure varies along the length of the evaporator 
tube as the liquid is heated to the boiling temperature in 
the non-boiling zone leading to two-phase flow of vapor 
and liquid mixture along the boiling zone. Pressure drop 
correlations for single phase liquid, would apply to non- 
boiling section. Pressure drop characteristics of the boiling 
zone would be based on the hydrodynamics of flow of two- 
phase mixture. Various factors contributing to pressure 
drop in the evaporator tube and correlations available for 
pressure drop are discussed in this section. 


2. 2. 2.1 Pressure Drop for ITon-Boiling Zone : 

In the non-boiling zone, frictional pressure drop over 
a differential section of the tube can. be determined using 


equation (2.1). 

tr . h!xa_e 

f 31 2 g c 


( 2 . 1 ) 


The effect of roughness characteristics of the tube on the 
value of friction factor can be determined from equations (2.2) 
and (2.3) for laminar and turbulent flow respectively. 

f = 16/U^ e 


( 2 . 2 ) 
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l/f 1/2 = -2 log [(e/3.7D) + (2.51/N Re -f l/2 ] (2.3) 

The above equations are modified to account for the effect 
of viscosity during non-iso thermal flow of the liquid to 
give equations (2.4) and (2.5) for laminar and turbulent 
regions respectively. 

Divide f by (i^/^) 0,38 (2.4) 

Divide f by (P-a/p-^) 0 * 1 ^ (2.5) 


2.2.22 Pressure drop for Boiling Zone ; 

Pressure drop will consist of three components - 
frictional drop, drop due to momentum change, and elevation 
pressure drop arising from the effect of gravitational field. 
Equation (2.6) gives the various components of the total 
pressure drop. 


(S£) 

^dZ' Total 


= ( 


dP 

dz 


/My 

l dZ ; i 


(M) 

l dZ ; e 


( 2 . 6 ) 


Several correlations [15,1-5] are available in literature 
for estimating the two-phase flow pressure drop, for specific 
conditions of system geometry, operating conditions and phase 
properties. In general the two models are available for 
the analysis of two phase flow: homogeneous model and slip 


model. 
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The basic assumptions[l6 ,21] of a homogeneous- model are; 

1. Equal linear velocity of vapor and liquid phase 

2. Thermodynamic equilibrium between phases. 

3. A single phase fraction factor applicable to two-phase 
flow [21 ]. 

Homogeneous model considers the two-phase as. a single phase 
possessing certain average properties of each phase. The 
difficulty with such a model lies in the lack of a suitable 
method for calculating the pressure gradient due to acceleration. 

The alternative slip flow model distinguishes between 
the vapor and liquid phases which are assumed to flow at 
different velocities but each at uniform velocity; [20]« 

This model differs from the homogeneous model in that there 
is slip between liquid and vapor. Momentum and energy balance 
equations are applied separately for each phase in deriving 
the model [Chapter 3]. 

Lockhart-Martinelli [l6 ] developed a procedure for 
calculating the frictional pressure gradient of an isothermal 
two-phase flow based on correlation of the data obtained from 
horizontal flow of air and various liquids (air -water, hydro- 
carbon-air) at atmospheric pressure. The expressions of two- 
phase frictional pressure drop were developed in terms of the 
contributions of the single phase pressure drop values. 

The results obtained by Lockhart-Martinelli for isothermal 
gas-liquid two-phase flow were modified by Martinelli-Nelson[22 ] 
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to apply in forced circulation boiling in vertical tubes. 

Another modification suggested by Levy [25 ] is the momentum 
exchange model based on the assumption that the momentum 
exchange between vapor and liquid phases tend to maintain 
the frictional and heat losses equal for the two phases. 

Levy’s correlation shows limited agreement for the range of 
variables considered by Muscettola [26]. This is attributed 
to the fact that this correlation does not take into account 
the dependence on mass flow rate, and diverges as the quality 
approaches one [26]. The various correiations developed 
for predicting two-phase flow pressure drop are based on 
data for gas-liquid flow and no generalized correlation is 
available for predicting pressure drop in two-phase vapor-liquid 
flow [ 15 ]. Dukler et al. [ 15 ] have systematically compared 
existing correlations [ 21 ] and found that Martinelli-Nelson 
correlation predicts frictional pressure drop within + 30 
per cent of measured values [22]. Martinelli-Nelson correlation 
was also found to be independent of mass flow rates, 

2.2.3 Heat Transfer Correlations ; 

2.2.31 Non-Boiling Zone ; 

In the non-boiling zone heat transfer to the liquid 
inside the tube occurs by forced convection. The vast amount 
of experimental data available in literature, have been 
satisfactorily correlated by Dittus-Boelter [17] equation (2.7), 
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valid for N Re > 10,000, and L/D ? 60 

Nu = 0.0278 Re 0 * 8 p r °’ 4 (^S) °* 14 (2.7) 

!J 'w 

With increase in liquor viscosity, the flow pattern inside 
the tube can be either in transition on laminar region. 

For fluids of high viscosity, Coulburn 117," 29] proposed 
equation (2.8) 

S t -Pr 2 / 3 = 0.023 Re“°* 2 (2.8) 


where 


q _ h. 

t C 


Stanton Number 


The equation (2.8) is recommended for hea,t transfer to 
fluids with temperature dependent properties and is applicable 
for 0.7 < Pr .< 160, N Re < 10,000, and L/D > 60. In 
equation (2.8) St is to be evaluated at the average fluid 


temperature , 


while N Re , 


Pr are determined using arithmatic 


moan of the wall and fluid temperatures. For a large viscosity 


of the process liquor the flow inside tube can become laminar. 
In streams flowing in viscous motion natural convection 
currents can have profound effect on the heat transfer. 

Studies of natural convection have shown that Grashoff number 


satisfactorily correlated a- heat transfer data for viscous 
flow. Grashoff Number introduces the added influence of 


natural convection currents contributing to heat transfer in 
forced convection in the viscous region. Eubank and Proctor [17') 
critically surveyed the available data for viscous flow and 
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arrived at the following empirical equation (2.9) . 

1^,. 0 .14* "* A 7r 1 V-^ 

(Hu) (--) = 1.75 [_(Gz) + 0.04 [ (D/1) (Gr)(N Re )] 0W5 J 

(2.9) 

Equation (2.9) can he used to determine the heat transfer 
coefficient for laminar flow region. 

2.2.32 Boiling Zone : 

Prediction of heat transfer coefficient for boiling 
of liquids in forced convection is very much of interest in 
practice, due to wide spread use of boiling reactors in 
nuclear power generation. But no fundamental theory is yet 
available to predict the boiling heat transfer coefficient 
for forced convective boiling. In order to achieve a meaningful 
correlation of heat transfer coefficient for two-phase vapor 
liquid. flow, as encountered in boiling, it is necessary to 
measure local heat fluxes and temperature differences. 

De.ngler and Addoms[i2], Guerrieri and Talty [31], Bennet [23] 
and Schrock and Grossman [32] have measured such coefficient 
for -the boiling of liquid in vertical tubes and annuli. 

In the forced convective region, Dengler and Addoms [12 J , 
and also Guerrieri and Talty correlate their results according 
to equation (2.10) by analogy to analysis of Lockhart and 
Martinelli. 
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E—- = A [ ] n (2.10) 

n IP *tt 

where h^p and h^p denote heat transfer coefficient for 
superficial liquid and two phase vapor liquid flow respectively, 
X. . - lockliart-Martinelli parameter [Appendix A] 


k-r - D 0(l-xl -C p P T 0,4 

0.023 i [-i PJ— ] 0 ’ 8 [ -^ 


( 2 . 11 ) 


Values of the coefficients in equation (2.10) reported by 
different authors, are as follows; 


Dengler and Addons ; A = 3.5 


n = 0.5 


Wright 


A = 2.72 n = 0.58 


Schrock and Grossman: A = 2.5 n = 0.75 

Collier and Pulling: A = 2.7 n = 0.7 

Departure of experimental results from these correlations in 
all cases was attributed to appreciable contribution by 
nucleate boiling towards heat transfer. These correlations 
do not take into consideration influence of nucleate boiling. 
Schrock and Grossman have proposed equation (2.12) which also 
includes the effect of nucleate boiling. 


= B [ 


0 ' \ 


+ A[ (^— ) n ] 


( 2 . 12 ) 


Q - he a.t flux 
G - total mass flux 

- heat of vaporization 
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B 


o 


"bubble Humber 


The following values of constants B, A, n have been reported 
by Schork-Gro ssrnan [ 32 ] and Wright [33]. 

Sehrock.-Gro ssman 5 B = 7.39x10^, A =1.3 x 10“ 4 , n=0.66 

Wright : B = 6.70 x 10 3 , A = 3.5 x 10~ 4 , n =0.66. 


Schrock find Grossman studied the following range of conditions. 
System; Water-steam in vertical stainless steel tube 
Tubes; 0.118 in I.D. x 15,30, 40 in long 
0.237 in I.D. x 15 and 30 in long 
0.432 in I.D. x 30 in long 
Pressure; 42-505 psia 
Quality; 5-57 per cent 


Schrock and Grossman obtained +30 per cent fit of the 
experimental data with the mentioned constants. 

It is expected that the contributions of the nucleate 
•boiling heat transfer and forced convective heat transfer 
vary with flow conditions. These variations have been 
investigated by Chen [ 34 ]. Chen obtained empirically the values 
of two dimensionless functions S and P, that allow for the 
variations in the nucleate boiling and forced convective 
components respectively. 


b-ipp k '^nucleate + ^forced convection 


(2.13) 
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This correlation agrees very well with wide range of 
experimental data [23»12].Chen uses Forsteir-Zuber ' s [ 35 ] 
nucleate boiling heat transfer correlation and the Dittus- 
Boelter forced convective heat transfer correlation. At the 
present time no generalized correlation is available for 
prediction of nucleate boiling heat transfer coefficient. 

2,3 Correlations Recommended for POE ; 

2.3.1 Pressure Drop Correlations ; 

Equation (2.1) can be used to predict pressure drop 
over non-boiling zone reliably. The equation (2.1) is 
independent of any system, as long as the liquid is Newtonian. 
Pressure drop correlations for boiling zone in PCE are not 
many. Martinelli-Nelson correlation has been tested for 
forced circulation boiling of water, and this predicts 
experimental data within + 30 per cent. This correlation 
can be recommended for boiling zone in PCE equation (A. 7 ). 

2.3.2 Heat Transfer Correlations ; 

Equations (2.7) to (2.9) can predict accurately heat 
transfer coefficient in non-boiling zone of FCE. Por 
predicting film coefficient in boiling zone in FCE, Schrock- 
G-rossman correlation is recommended. This correlation has 

been developed for forced circulation vaporization of water 

$ 

in stainless steel tubes. This correlation takes into account 
both the mechanisms of heat transfer in the boiling zone. 
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The range of heat flux and pressure considered by Schrock- 
Grossman agrees with FCE conditions. This correlation is 
expected to predict film coefficient for boiling zone well. 

2.3.3 Recommendations for Kraft Black Liquor System : 

Present work considers kraft black liquor as the 
process fluid in the POE unit. Black liquor was selected 
since the correlations for various engineering properties are 
readily available from previous work [36,37] .The various 
correlations discussed earlier have been developed for systems 
different from black liquor. Comments on the applicability 
of these correlations for the black liquor system are given 
in this section. 

Pressure drop for non-boiling zone can be reliably 
predicted from equation (2.1). This equation is dependent 
only on density, friction factor Reynold's number. Film 
coefficient for black liquor in non-boiling zone can be 
predicted from equations (2.7) to (2.9). Harvin [17] has 
worked extensively with black liquor system and verified the 
validity of the equations (2.7) to (2.9) for predicting film 
coefficient in horizontal tubes for laminar transition and 
turbulent flow conditions for the concentration range of BO to' 
60 per cent solids. 

Martinelli-Nelson correlation, based on experimental 
data for forced convective vaporization of water and Lockhart- 
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Martinelli parameter for frictional pressure drop, gave 
satisfactorily prediction for the pressure range of 120 to 
20000 kl\f/n 2 [22] and quality 4 to 95 per cent. Because of 
the applicability of Martinelli -Kelson over a wide range of 
pressure and quality it is recommended for black liquor 
system even tho ,gh the effect of viscosity is not included in 
Martinelli-Kelson correlation. It is assumed that Martinelli- 
Kolson equation will give a satisfactory estimate of pressure 
drop, over a moderate range of values. 

Prediction of film coefficient in boiling zone also 
raises the question of applicability of Sohro ck- Grossman corre- 
lation for black liquor system. Even though the correlation 
considers both nucleate and forced convective mechanisms, for 
water, the effect of viscosity is not well defined. A complex 
system of like black liquor is not convenient in fundamental 
studios of 2-phase flow heat transfer/fluid flow mechanisms. 

In the present work the viscosity range is 4 to 30 cp and 
the liquid is Newtonian. Schrock -Grossman correlation will 
overestimate the heat transfer coefficient by probably 30 
per cent. ‘This limitation must be borne in mind while 
interpreting the results obtained by the simulation model 
of this study for the black liquor system. 



CHAPTER 5 


MATHEMATICAL MODEL OP FORCED CIRCULATION EVAPORATOR 

Simulation requires simultaneous mathematical description 
of the hydrodynamic and heat transfer phenomena occurring 
along the length of the evaporator tube. The present study 
deals mainly with steady state simulation. Mathematical 
model for the evaporator is developed from basic energy, 
continuity and momentum balance equations and correlations 
available in current literature for heat transfer in non- 
boiling and boiling zones of the evaporator tube. Simulation 
of an evaporator unit consists in predicting longitudinal 
profiles of liquor temperature, concentration, heat flux, 
overall heat transfer coefficients and evaporation rate. 

To obtain overall performance a differential approach is 
adopted by considering the tube length to consist of a 
series of differential segments (dz) . Momentum, energy 
balance, continuity equations together with various corre- 
lations are applied for a differential length (dz) to obtain 
mathematical expressions which predict changes in liquor 
conditions over the differential length (dz) by an iterative 
procedure. The desired longitudinal profiles are then 
obtained by a sequential traverse of the non-boiling and 
boiling zones of the entire tube length. The model is developed 
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for ail evaporator with vertical calendria and the modifications 
necessary for a horizontal calendria are discussed in Section 3.4 

3.1 ICE with Vertical Calendria ; 

3.1.1 Eon-Boiling Zone ; 

Energy, momentum and continuity equations are applied 
over a differential length (dz) for the one dimensional 
analysis of the single phase flow in the non-boiling section 
(Eig. 3.1a). The following assumptions are used in the 
analy sis. 

1. Radial variations in parameters like pressure , density , 
temperature. Velocity and other transport properties 
are considered to be small* 

2. Steam is at saturation temperature, 

3. Heat losses are negligible, 

3.1.2 Momentum Balance : 

Equation (3.l) gives the momentum balance for the 
differential section. 

= Rate of change of momentum of fluid (3.1) 

^p ^p ^ 

Momentum input = — up 

s c 

( ^p d $ p ) ( Up + dUp ) A 

Momentum output = (%, + du™) 

°c 

Change in momentum = J; [ ( /p+d /p) (up+dup) 2 - ^pUp] 

(3.2) 
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&.2TTQ dZ 


FI GUKE 3.1a; DIFFERENTIAL BE GHENT IN NON-BOILING ZONE 



FIGURE 3.1b; DIFFERENTIA SEGMENT IN BOILING ZONE 
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Equation (3.3) represents material balance 
Afp Uj, = Constant 
1 (-A- jfp n.p) = 0 

Jp du-p = ~UpdJ(p (3.3) 

Combining equations (3.2) and (3.3) gives change in momentum; 

g c (3.4) 

Equation (3.5) gives net forces acting upward on the fluid 
element and include contributions from pressure, hydrostatic 
head and shear .forces. 


A.P ~ A*(P + dP) - JL-A-dz- f - t •ifD.^dz (3.5) 

Equation (3.6) relates frictional pressure drop and shear 
stress at the wall. 



Net upward force on the fluid element is obtained by combining 
equations(3.5) and (3.6) to give (3.7) 


-A*dP - 



A»dz 


(3.7) 


Substitution of the above expressions for the various terms 
of equation (3.1) results in equation (3.8) fbr the prediction 
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of pressure gradient along the length for the non-boiling 


zone . 


dP 

dz 


i~s _ (Ik) . nil du ? 

, F g c dz ' g c dz 


(3.8) 


3.1.3 Total Energy Balance ; 

The general form of total energy balance equation is 
as follows: 

d (Enthalpy) + d (Kinetic Energy). + d (Potential Energy) 

= dQ ' (3.9) 

The calculation of each term in equation (3.9) for the 
differential section is shown below. 

d (Enthalpy) = Enthalpy of liquid leaving - Enthalpy of 

liquid entering 

- ; _ = W T t(C p + dCp) (T+dT - T q ) - (^(T-T^)] 


d( Kinetic Energy) 


dT + dCj, (T-T q )] 

1 „ r (u P + au P )2 3 . 

2 W T t g c J g 0 J 

- Up du p 

V r 1 -] 


] 


Sc J 


d (Potential Energy) = • 'dz 


Replacing individual terms in equation (3.9) and assuming a 
constant value for specific heat (Cp) gives equation (3.10) 



for predicting the change in temperature of liquor (dT), 


w t[°f dT + ~i-J 


Up * 


+ 


-g_ 




dz] = Q_^ 71 D^ dz 


dT 


■'F 


- Q. it D. dz 
[— X 


u-p du p 

- fS I®] 


% 


T 


So J 


Sr J 


( 3 - 10 ) 


3.1*4 Boiling Zone : 

The following simplifying assumptions are made for 
developing necessary equations applicable to the boiling zone, 


1. Plow is one dimensional 

2. Negligible radial variations of velocity, pressure, 
temperature and quality. 

3. Equilibrium exists between vapor and liquid in the 
boiling zone. 

4. Density, viscosity, and other transport properties 
are taken as weighted average of the mixture. 


3.1.5 Momentum Balance ; 

Application of momentum balance equation (3.1) to a 
differential section dz in the boiling zone will give (3.11) 
for the total pressure gradient (Figure 3.1b) 


,dPx _ /dP\ _ /dP\ _ ,dl\ 
'dz t dz'f ” uz'a ” Mzm 


(3.11) 


Nlz't 



= Total pressure gradient 

= Frictional pressure gradient obtained by 
Martinelli-Eel son correlation [22]. 
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/dPN _ 

IktT/ e ~ Pressure gradient due to elevation change 


l dz ; e 


(dP) 

% a 


Accelerative pressure gradient (discussed in 
Appendix A) 




d_ 

dz 


[ i 


1-x) 2 x 2 

J p (l-ra) . j > v a 


] 


( 3 . 12 ) 


Density and viscosity of the two phase' mixture are taken 
as the weighted average of the values of the separate phases 
using quality (x) as represented by. equation (3.13) 


i/f 

i/f 


1 - x 
—rr — 

J F 

1 - x 

Up 


x 

r 

J v 


p. 


(3.13a) 

(3.13h) 


v 


3.1.6 Total Energy Balarice : 

Enthalpy of mixture leaving dz: 

W T (l-x-dx) Cj,(T+dT^,T 0 ) + W t (x-kLx) C T ( T+dT-T Q ) 

+ ¥ T (x+dx)A v 

Enthalpy of mixture entering dz: 

W T (l-x)» (T-T q ) + ¥ T C v x (T-T q ) + X v 

d (Enthalpy) = ¥ T * Cptl-x) dT - (T-T Q )*dx + ¥ T „C v x.dT 


+ V/ T .C v (T-T Q )dx + ¥ T *dx A 


V 
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a 


(Kinetic Energy) = A. d[ + 

J Y * 


Substituting the individual terms in the total energy 
equation gives equation (3.14) 


dT 


1 


Cp (l-x) + C y x C 
-dx/S 


" Q i * D i dz 


' ¥ 


T 


+ C F (T-T 0 )dx - C v (T-T Q )dx 


r rjt 


d (y^ + 




) - dz ] 


v 2g c ^ a j>2 (1 _ a) 2 g c J 

(3.14) 

Equation (3.14) can be rearranged to give equation (3.15) 
for predicting the change in quality ’dx' , by choosing 1 as 
the reference temperature. 


dx 


A 


1 -Q.it D. dz 

1 [ 1 3- 


G; 


-X- d r— Si + jlzSil 

2 g 0 J L ? a a 2 f V (1 _ a) 


V 


w, 




T 


] 


- §- T dz 

S C J 


- dT [C v x + Cp(l-x)]] 


J? 


(3.15) 


3.2 Mass Balances : 


3.2.1 Solute Balance in Boiling Zone : 

Solute balance over the differential section gives 
equation (3.16) for predicting changes in liquor concentration. 
W T (l-x) S = W T (l-x-dx) (S + dS) 

ds (3 - 16) 
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3.2.2 General Mass Balance for Evaporator with Recycle ; 

Figure 3.2 gives different variables associated with 
the streams for an evaporator arranged for recycle of a 
portion of the product to accomplish the necessary rise 
in liquor concentration. 

Total material balance : 

\ = W 1 - W 2 (3.17) 

Total Solid balances 


w i s i 


^ 2^2 


Recycle ratio: 

¥ 

R _ i 
R ~ ¥ 2 


(3.18) 


(3.19) 


Total food to evaporator: 


¥ 


T 


W 1 + W R = w i + aw 2 
W 1 S 1 

h + ~iy R 

s 

¥ 1 (l+"| i '• ' R) 


( 3 . 20 ) 


Inlet liquor concentration to the evaporat&r tube. 


S, 


¥ 1 S 1 + ¥ R S 2 
= ¥ x + ¥ r 


¥- l S 1 (l+R) 


¥ 


T 


(3.21) 


3.3 Material and Energy Balance for Flash Separator s 

Liquid issuing from tubes is flashed in vapor space 
and achieves an additional increase in liquor concentration. 
Vapor space is invariably maintained under vacuum. 
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and connected to an intermediate effect of a multiple 
effect evaporation plant. A large chamber is used for 
flashing to reduce the velocity of vapor stream and to 
reduce entrainment of liquor droplets. 

Temperature of liquor leaving flash tank (3? 2 ) is 
taken as reference temperature for energy balance equations 
of the separator. 


Equations (3.22) and (3.23) represent inlet (H^) and 
exit enthalpies of the vapour -liquid mixture. 

H i = w i cpVh) + W v °v (W + K (3.22) 

Ii 0 = \ \" v + ¥ v C v (T 2 -T 2 ) + ¥» C” (T 2 -T 2 ) (6.23) 

\ = % + W » 

Enthalpy balance gives equation (3.24) 


Of (Tt-T 9 ) + (T-^) + v (3.24) 


w 1 

W L v ' i 'l _J '2 


Total material- balance around flash separation 


VL 


¥ * + ¥ • 
V Jj 




Substituting for ¥ v in enthalpy balance gives equation (3.25) 
for the rate of flash liquor , , 

w;[ X 1 ; -\' v - c;(-J 1 -i 2 )]t- w-[ T - cpyrp] 


w i’ 


X 


// 


(3.25) 


Liquor concentration after flashing is calculated from 
equation (3.18) 
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3.4 Modif ications for Evaporator with Horizontal Qalendria : 


Evaporations with, horizontal calendria, have a horizontal 
heating element. To predict the necessary profiles, elevation 
terms are deleted from equations (3.8), (3.10), (3.1l) and 
(3.15). Equations (3.8) and (3.10) for the non boiling zone 


are arranged to give equations (3.26) and (3.27) 


^dP^ U E % du E 

\Iz ; ~ Mz ' g 


(3.26) 


dT 


Op 


[' 


Qi 71 D i dz 


1 F 


T 


Up dup 


■] 


(3.27) 


For the boiling zone, equation (3.15) is modified to equation 
(3.28) for predicting change in quality (dx). 


-CL % D. dz 
dx = [. 1 1 


g; 


A 


v 


¥ 


T 


2g c J 


1 a ( - z 


(1-x) 


) 


n '« 2 

-dT *[0j,(l-x) + C v x]] (3.28) 


Simulation computations are carried out in a stepwise manner 
assuming a series of differential sections of the evaporator 
tube, and using the criteria of maxima in liquor temperature 
profile for the transition from non-boiling to boiling zone[lo ] 

3.5 Engineering Properties of Process Liquor : 

Computations require correlations for engineering 
properties process liquors. Necessary engineering data and 
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correlating are available for kraft black liquor system 
from earlier work £36,37]. 

(a) Density (g/cc) : 

% L = 1.0182 + 0.779 • S - (0.000219 + 0.0003 • S(*T 

(3.29) 

(b) Thermal conductivity (kcal/hr.m °C) 

k L = 0.53354269 - 0.3425-S + (0.0011355 + 0.000209-S) -T 

(3.30) 

(c) Specific Heat(kcol/kg °C) 

C L = 0.993 ~ 1.25*S + (0.0000506 + 0.0104*S)*T 

(3.31) 


(d) Do generalized correlation is available for viscosity 
of black liquor in the concentration range 35-60 per cent. 
Table 3.1 gives viscosity data for bamboo kraft black liquor 
over a range of concentration and temperature values. from 
the available experimental data a correlation of the form 



(3.32) 


has boon d eveloped by polynomial regression analysis, for 
each reported concentration of black liquor. Viscosity for 
a concentration of black liquor, intermediate between two 
successive concentrations of Table 3.1 is estimated by 
linear interpolation. 



36 


TABLE 3.1 

VISCOSITY BAMBOO ZBAFT BLACK LIQUOR 


Percent 

solid 



Temperature , 

°C 



50 

60 

70 

80 

90 

100 

iio 

35 

10.23 

7.71 

. 5.80 

. .-.4.37 

3.29 

2,48 

1.88 

40 

21.58 

14.64 

10.45 

7.81 

6.06 

4.87 

4.03 

45 

40.26 

25.66 

16.77 

11.22 

7.66 

5.33 

3.78 

50 

68.27 

41.52 

27.23 

14.06 

14.11. 

10.98 

8.91 

55 

173.97 

99.00 

61.63 

41.58 

29.98 

22.93 

18.45 


TAPLE 5.2 


PRO PERTIES QE LIQUID WATER 


Temperature 

°c 

Viscosity 

cp 

Thermal Conducti- 
vity ? kcal/h m °C 

"Density 

gn/cc 

38 

0.682 

0.539 

0.994 

49 

0.559 

0.552 

0.989 

60 

0.470 

0.563 

0.984 

71 

0.401 

0.572 

0.978 

82 

0.347 

0.578 

0.971 , 

93 

0.305 

0.584 

0.964 

104 

0.270 

0.587 

0.954 
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(g) Boiling Point Elevation (°C) 

BPE = (-2.57709 + 30.7781-S + 1.49186*S 3 )/1.8 (3.33) 

(f) Vapor Pressure - Temperature Relationship for Water; 

Vapor pressure and temperature can be correlated by 
Antoine equations (3.14) and (3.15) 

l0 « p = A - dk+c 0.34) 

T = £ C + BPR (3.35) 

A-logP 

(g) Properties of liquid water as a function of temperature 
has been generated from Table 5.2 by linear interpolation. 

3.6 Void Fraction Correlation. ; 

In the boiling zone correlation for predicting void 
fraction is available from literature. Anthony etal have obtained 
correlation by regression analysis of large volume of data [1-9]. 
The correlation is given by equation (3.36). 

a = ©xp[-0. 25522 - 0.10583 • ln(X) - 0.0289 • (ln(X)) 2 

-0.00884 • (ln(X) ) 2 ] (3.36) 

a - fraction of volume occupied by vapor in a mixture of 
vapor and liquid. 

X - Lockhart Mnrtinelli Parameter. 
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3.7 Co mputational Scheme ; 

Simulation computations are carried out in a stepwise 
manner assuming a series of differential sections of the 
evaporation tube and using a criteria of maxima in liquor 
temperature profile for transition from non-boiling to 
boiling zone. 

3.7*1 Computational Algorithm for Non— Boiling Zone ; 

For a weak black liquor feed below its saturation 
temperature, flow chart of computation for the non-boiling 
section is given in Figure 3 . 3 . Flow rate, concentration 
temperature, pressure, recycle ratio, steam chest temperature, 
and tube geometry are fixed for starting calculations. A 
differential length (dz) of tube is considered to determine 
the heat transfer rates. Temperature differential dT is 
assumed first. Liquor properties at mean temperature, T+dT/2 
is evaluated. The change in pressure in this section is 
calculated from the frictional pressure drop correlation 
equations (3.8), '(2.1 ), (2.3), (2.5), (3*26)- 

The overall heat transfer coefficient for the differential 
section is calculated from the individual film coefficients 
for liquid and steam sides and resistances of the wall and 
scalo deposits. Heat flux is computed from heat transfer 
rate equation (Q = U.AT ) . Equation (3.10 ), (3. 27) are 
used to determine differential change in temperature dT and 
compared with initial assumption. If the calculated value 




FIGURE 3.3: COMPUTATIONAL ALGORITHM FOR NON-BOILING SECTION 
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and assumed value differ by more than a preassigned limit, 
calculation is repeated with new value of dT. Otherwise 
calculated liquor temperature is compared with saturation 
temperature, corresponding to pressure at the exit of the 
differential section. If the liquid temperature happens to 
bo below the saturation temperature, calculation is continued 
for next section with a new assumed value of dT, This 
process is repeated until the liquor temperature equals the 
saturation temperature of liquor. This point corresponds to 
the end of non-boiling regime, and subsequent calculations 
are carried out for boiling section. 

3*7.2 Computational Algorithm for Boiling Zone. : 

* . 

The whole tube length is divided into a fixed number 
of segments. Computational algorithm for estimating changes 
in the process liquor conditions in the boiling section is 
given in Figure 3.4. A differential section of length 
(dz = 0.1 m) is considered in the boiling section. Calculations 
over the differential section give the change in liquor to 
temperature, pressure, pressure gradient, concentration, 
quality, void fraction heat flux and overall heat transfer 
coefficient. Value for the change in dT is assumed for the 
differential section. Pressure gradient over the differential 
section is calculated using equations ( 3.11 ) to ( 3.13 ) an( ^ 
correlation given in Appendix A. Total pressure drop over the 
differential section is obtained by multiplying pressure 





FIGURE 3.4: COMPUTATION ALGORITHM FOR BOILING SECTION 










43 


gradient Toy dz. The temperature at the end of th e exit 
of the differential section is calculated from vapor-liquid 
equilibrium for water (Appendix — c ), corresponding to the 
pressure at the exit of the differential section. To this 
temperature boiling point rise ; is included. - Overall heat 
transfer coefficient is calculated from individual film 
coefficients, and contributions of the tube wall and dirt 
resistances. Heat flux Q is calculated from Q=U*A T and 
equation (3.1^is used to calculate change in quality. It 
is assumed that vapor and liquid remain in equilibrium, 
and calculated dl is compared with assumed dT for convergence. 
When temperature difference dT has been determined, change 
in quality (dx) , and liquor concentrations (dS) are determined 
from equations. ( 3 . 15 ), ( 3 . 16 ) A counter "’Z’ keeps track of 
number of segments calculated and calculation is stopped when ' K* 
reaches 'NSEC total number of segments in a pass. Similarly 
calculation is repeated for next pass, and continued until 
calculation for the- all passes are completed. 

3.7.3 Computational Algorithm for Plash Calculation: 

In FCE major portion of evaporation takes place in 
the vapor separator and. the change in liquor concentration 
is computed according to Figure 3 .5 . Initially outlet liquor 
concentration Sg is assumed. From equation (3.18) flow rate 
of concentrated liquor is calculated. For a given recycle 













45 


ratio, equations (3.19) to (3.21) are used to calculate 
total flow rate and inlet liquor concentration to evaporator, 

feed temperature and feed inlet pressure are known. With, the 
given inlet conditions algorithms in Figure 3.3 and Figure 3.4 
are used to determine liquor conditions and the end of tube, 
liquor concentration after flashing is calculated using 
equations (3.25) and (3.18). Outlet liquor flow rate (W» ) is 
calculated from equation (3.25). Liquor flow rate to flash 
separator (WjO , and inlet liquor concentration are known. ■ So 
outlet liquor concentration can be calculated from equation 


(3.18). Now the assumed value of S 2 and calculated value of 


S 2 , S 2cal arG compared. 


The calculated value is compared 


with the initial assumption (S 2 ) for convergence. 




CHAPTER 4 


RESULTS AMD DISCUSSIONS 

Computational schemes outlined in Figures 3.3 to 
3.5 are used to simulate the evaporation process-. The 
model is used for the following three cases: . 

1. Simulation of Badger's experimental data [10] 

2. Simulation of forced circulation evaporation 

of kraft black liquor in a single: vertical 
,tube. - ■ 

3. Simulation of forced circulation evaporation 

of kraft black liquor in a single vertical/ 

' . , . , ■„> . ... . ., . , - ■ 1 . : ; ■ '■ . ■ ■ ' '■ . . .■ • ' . ■ . . ' 

horizontal tube with recirculation of liquor. 

The results of the simulation studies for the above 
cases are discussed in the following sections. 

4.1 Simulation of Badger’s Experimental Lata [ 10] . 

Badger and Brooks [10] have reported experimental 
data for boiling heat transfer obtained from a long tube 
vertical unit equipped with a single tube (L=6.1 m, 44.5/ 
50.8 mm dia.) using water as the system. The above data 
include boiling length, evaporation rate, overall heat 
transfer coefficient and total heat transferred in boiling 
zone. Seventeen sets of data were drawn from Badger and 
Brooks' study and used in this investigation and are given 
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TABLE 4.2b 

QOjjMi^LSOIi OF BADGER'S DATA WITH SIMULATED RESULT 


SI. 

Eo, 

Obs. 

E ( ks:/h) 
Cal. 

Percent 

Dev. 

Obs 

A b Q b (kW) 
Cal* 

Percent 

Dev. 

1 

59.69 

80.10 

25.48 

, 37.68 

50.32 

25.12 

2 

54.79 

67.75 , 

19.12 

34.74 

' 42.53 

18.31 

3 

56.70 

76.15 

25.54 

36.17 

47.16 

23.30 

4 

73.48 

100.43 

26.88 

47.54 

62.96 

24.49 

5 

63.50 

73.39 

13.47 

39.55 

47.66 

17.01 

6 

66.86 

79.70 

16.11 

41.84 

48.67 

14.03 

7 

29.30 

30.20 

■ 2.98 

18.00 

18.94 

4.96 

' 8 

41.00 

52.99 , 

22.62 

. ... • S ¥ 

25.65 

33.76 

24.02 

■9 

25.31 

21.18 

-19.49 

15.34 

! f * ■ 

13.24 

-15.86 

10 

81.46 

166.94 

51.20 

50.92 

101.45 

49.80 

11 

20*77 

23.21 

10.51 

13.69 

14.37 

4.73 

12 

38.96 

47.98 

18.79 

24.68 

29.79 

17.15 

13 

38.62 

26.23 

-47.23 

17.50 

15.52 

-14.98 

14 

43.90 

49.73 

11.72 

27.51 

30.07 

8.51 

15 

34.33 

37.54 

8.55 

21.00 

21.47 

2.18 

16 

34.20 

37.08 

8.65 

20.74 

21.43 

3.21 

17 

31.29 

41.46 

24.52 

18.56 

17.35 

-6.97 


1,1; . • : : FUH 

L 1 *- »?ARV 

~~Za 5835a 

Am* No. £% ** 
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in Table 4.1. Numbers in parentheses correspond to original 
serial number in Badger and Brooks' experimental data. The 
length of differential se'ction (dz) was taken as 15 cm for 
the computation. The feed was below saturation temperature 
in all the cases giving both non-boiling and boiling zones. 

The length of non-boiling zone was predicted using equations 
(2.5) to (2.5), (2.7) to ( 2. 9) , (3.3) and (3.10) by stepwise 
computation till the maximum in temperature profile was 
obtained. Average values of film heat transfer coefficient 
and overall heat transfer coefficient were calculated from 

h inb = 2 h inb * ^ T i // 2 T i and U nb = 2~ i >AT i / ^ T i* 

Subsequent calculations based on equations (2.6), (2.10), 

( 3 . 10 ), ( 3 . 12 ), (3.15) gave liquid film, overall heat transfer 
coefficients, evaporation rate and total heat transferred 
in boiling zone. Table 4.2a and 4.2b give a comparison of the 
above results with Badger and Brooks' experimental observations. 
The results show very good agreement of predicted length of 
non-boiling zone with the observed values with an average 5 
per cent deviation. Predicted values of overall boiling heat 
transfer coefficient, evaporation rate and total heat trans- 
ferred in boiling section show average deviations of 13, 13 
and 12 per cent respectively. Earlier simulation study by 
Agarwal [9] has also considered some of the results of Badger's 
and Brooks. The average deviations obtained in this study 
and values reported by Agarwal [9] are given below. 
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Deviation, Per cent 


Variable 

This Study 

Agarwal 

L b (m) 

5 

15 

U b (kW/m 2 K) 

13 

13 

E(kg/h) 

13 

13 

A b Q b 

12 

17 


The above values show a better agreement of the 
prediction of this study with Badger’s observations. Liquid 
film heat transfer coefficient in non-boiling zone is predicted 
in this work by Eubank and Proctor [ 30 ] equation for the 
viscous region which includes Grashoff number to account for 
the effect of natural convection. Agarwal has used an 
empirical correlation [40 ] for estimating the overall heat 
transfer coefficient in non-boiling zone which leads to reported 
15 per cent average deviation for the length of non-boiling 
zone, fable 4.3 shows comparison of Badger's result with 
Agarwal ' s result. 

Heat transfer coefficients for boiling in Agarwal' s[ 9 ] 
model are based on Coulson-McNelly and Penman-Tait [11,14] 
correlations for the nucleate and annular film boiling regions 
respectively. Schrock-Grossman correlation for forced 
convection boiling is used in this study to predict film heat 
transfer coefficient for the entire boiling zone. The two 
models predict Badger' s data equally well as observed from 
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the average deviations for the value of boiling length, 
overall heat transfer coefficient, evaporation rate and total 
heat transferred in boiling zone. 

A summary of the simulation results of this study 
using Badger* s experimental data is given in fable 4.4. Table 
4.4 shows temperature and pressure of t he fluid at the inlet, 
end of non-boiling zone and exit of the tube, quality and void 
fraction of vapor liquid mixture at exit and liquid film and 
overall heat transfer coefficients in boiling zone. A typical 
sot of profiles for temperature, pressure, quality, film and 
overall heat transfer coefficients are given in figure 4.1. 
for data set 11. Temperature profile shows a maximum of 70.5°C 
at 1.65 m indicating completion of non-boiling zone, this is 
accompanied by a rapid reduction in pressure caused by decrease 
in hydrostatic head. The values of film and overall heat 
transfer coefficient are nearly constant at 1.05 and 0.90 
respectively. The start of boiling zone is evident from 
figure 4.1, from the profiles of quality, and heat transfer 
coefficients as depicted by the graphs in figure 4.1. These 
profiles are similar to the results obtained by Agarwal. 

Thus both the models satisfactorily predict Badgers and Brooks 
o xp e r im o n t al data. 
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4 • 2 J?jL.jfQ- r PQd Circulation Evaporation in a Single 

Vortical Tube ; 

The effect of various operational parameters like 
feed temperature, feed liquor concentration, feed rate on 
performance of a single tube of 6 m long vertical evaporator 
is studied using the simulation programme. Even though in 
commercial practice evaporator units consists of 200-500 tubes 
(single pass) in the calendria, in this work, a single tube 

is considered for simulation purposes. The single tube is 

. , , representing 

assumed to do a typical one^fluid flow and heat transfer 

mechanisms encountered in a full scale unit* 

It is assumed that liquor is distributed equally in all the 

tubes, and that the various process conditions are identical 

for all the tubes. A dirt scale resistance = 0.45 m 2 K/kW is 

assumed. 

In this study, black liquor is taken as the process 
liquor, because the necessary correlations for engineering 
properties for various commercial samples |a«.a? ] are readily 
available, further forced circulation evaporator units are 
increasingly used for the viscous concentration range of 
55-65 per cent. 

4.2.1 Effec t o f fe ed Liquor Concentration on Performance ; 

feed liquor concentration of bamboo kraft black liquor 
is varied from 35 to 65 per cert, while keeping the feed 
temperature and rate constant at 100°C and 10 kg/s/tube. The 
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latter corresponds to a tube velocity 3.5 - 4.0 m/s. 

Simulation results of Table 4.5 show that the overall heat 
transfer coefficient falls with increase in feed liquor 
concentration which contributes to a higher viscosity. viscosity of 
liquor increases 2.5 cp for 35 per cent to 48 cp for at 60 
per cent concentration (at 100°C) with a decrease in N Re from 
10x10^ to 0.5 x 10^. Figure 4.2 shows the profiles of non- 
boiling and boiling overall heat transfer coefficients with 
liquor concentration as the parameter,. Figure 4.5 shows average 
film heat transfer coefficients for non-boiling and boiling 
zones and overall heat transfer coefficient decrease with 
increase in viscosity of feed. Pressure profiles along the 
tube are depicted in Figure 4.4 with feed concentration a s 
parameter show that increasing feed viscosity pressure at the 
tube end is lower. Effect of feed liquor concentration, feed 
liquor temperature and feed rate on total pressure drop is 
given in Figure 4.5. Total pressure drop as depicted in 
Figure 4.5 increases with the viscosity of the feed. The 
results in Table 4.5 also show that change in liquor concen- 
tration is very small (less then 1.0 percent) and this 
necessiates the use of recycle concept to achieve the desired 
concentration levels. 

4.2.2 Effect of Feed Temperature on Performance ; 

The effect of temperature pjf bamboo kraft black liquor 
(45 per cent) was studied at 90 to 110 C and the predicted 
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results are summarised in Table 4.6. The results show an 
increase in overall heat transfer coefficient and decrease 
in pressure drop with feed temperature. Figure 4.6 shows 
profiles of tube side film coefficient for various feed 
temperatures.. There is a discontinuity in the profile during 
transition from non-boiling to boiling zone. Film coefficient 
for boiling zone slightly falls below non-boiling zone and 
then rises. Film heat transfer coefficient for boiling- zone 
is always higher than film heat transfer coefficient in non- 
boiling zone. This slight discontinuity can be due to failure 
of Schrock-Grossman correlation to predict film heat transfer 
coefficient during transition from non-boiling to boiling 
region. Schrock-Grossman correlation is developed on the basis 
of data for fully d evel oped boiling zone. 

Figure 4.7 shows an increase in average film heat 
transfer coefficients for both non-boiling and boiling zones. 
This is due to the fact that feed liquor concentration remaining 
same feed liquor viscosity falls with rise in feed temperature 
and consequently inlet Reynold's number increases. Overall heat 
transfer coefficients also improves due to increase in the 
individual heat transfer coefficients. 

Table 4.6 shows that liquor concentration at the exit of 
the tube increases with increase in feed liquor temperature. 
Figure 4.8 shows the concentration profiles with feed tempera- 
ture as parameter. Increase in liquor concentration at the 
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exit is due to the fact that with increase in feed 
temperature, boiling length increases and more vaporization 
takes place inside the tube. For a feed temperature of 90° C 
there is no boiling section and liquor concentration remains 
unchanged. For a feed temperature of 110°C the boiling- 
length is 3.4 m and concentration rises by 0.3 per cent. 

Figure 4.9 shows pressure profile along tube length indicating 
that with increase in feed temperature boiling length increases a r 
Figure 4.5 shows a marked fall in pressure drop with increase 
in feed liquor temperature. This is due to fall in viscosity 
of the feed liquor with increase in feed temperatures. 

4.2.3 Effect of Feed Liquor Rate on Performance ; 

Feed liquor rate of black liquor was varied from 5 to 10 
kg/s/tube, while feed liquor concentration was kept fixed at 
40 per cent and feed temperature at 100° G. 

Table 4.7 summarises the predicted results. Figure 4.10 
shows profiles of tube side film coefficient for various feed 
rates. Slight discontinuity in the film coefficient has been 
explained in Sec. 4. 2. 2. Figure 4.11 depicts effect of feed 
rate on average heat transfer coefficients. Table 4.7 shows 
10 per cent increase in heat transfer coefficients and 15 per 
cent increase in pressure drop for a two fold increase in 
feed rate. 
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4-3 Ilo iffed Oiroulation Evaporation with Recycle ; 

Concentration rise per pass in a PCS unit without 
recycle is very small (less than 1.0 per cent), hence to 
achieve higher concentration of the product liquor, part of 
the product is recycled. A commercial FCE unit normally 2—6 
passes and .100-200 tubes/pass. An evaporator consisting of 
three passes and one tube/pass is considered. An evaporator 
consisting of three passes and one tube per pass is considered 
in this work. Effect of feed liquor temperature, feed liquor 
concentration and recjrcle ratio on performance are studied 
using simulation programme. Results for vertical and horizontal 
orientations of the calendria are compared. Kraft black liquor 
is taken as the process liquor. Vacuum in the vapor space is 
assumed to bo constant at 35 kii/m 2 . 

4.3.1 Effe ct of Peed 'te mperature on Performance : 

Peed temperature is an important process variable in 
evaporator operation and determines effective utilization of 
the available heat transfer surface for evaporation duty. Peed 
temperature of inlet liquor is varied (90, 100, 110° C) keeping 
other parameters constant to study change in performance. 

The values of feed concentration, recycle rate, and separator 
pressure are kept constant. Tables 4.3 and 4.9 summarises 
the results for POE with vertical and horizontal calendria 
respectively. The results show a marked increase in product 
concentration with increase infeed temperature. Peed- liquor 
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with ail ini bial concentration of 45 per;. cent gives final 
product concentration range of 49 to 56 per cent, when 
temperrature of. feed liquor changes from 90 to 110°C. for 
evaporator with horizontal calendria, most of the evaporation 
takes place in the separator and no boiling occurs inside 
calendria. figure 4.13 shows -.a slight fall in overall heat 
transfer coefficient with increase in feed liquor temperature. 
This can be explained by ...the fact that, with increase in feed 
temperature concentration of mixed feed (fresh feed plus 
recycle liquor) increases and the resultant increase in 
viscosity lowers and in turn causes a decrease in overall 

heat transfer coefficient, figure 4.12 shows an increase in 
overall heat transfer coefficient with increase in feed liquor 
temperature. The latter increases the length of the boiling 
zone and gives a higher overall heat transfer coefficient." 

4.3.2 E ffe c t of fee d L iquor Concentration ; 

feed liquor concentration is varied from 35 to 55 per 
cent, while keeping other parameters, feed temperature, recycle 
ratio, vacuum in the separator constant. Tables 4.10 and 4.11 
summarises the simulation results for evaporator with vertical 
and horizontal calendria respectively. It is seen from 
the Tables 4.10 and 4.11 that overall heat transfer coefficient 
falls with increased concentration which raises the viscosity. 
Viscosity of liquor at 35 per cent concentration, increases 
from 4.8 cp, 100°C to 52 cp at '55 per cent concentration and 
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100 G. lor evaporator with vertical ealendria boiling 
occurs in the tubes and the overall heat transfer coefficient 
is higher for all concentration than non-boiling, heat . transfer 
coefficient for evaporator with horizontal ealendria. 
figure* 4.14 and figure 4.15 show the overall performance 
based on Table 4.10 and 4.11. Total pressure drop increases 
with feed concentration and higher for evaporator with 
vortical ealendria, due to the contribution of elevations 
head loss. 

4 .3. 3 Affect of Ho cycle Ratio : 

Tho rise in concentration of the ; liquor obtainable from 
a i'CTl unit with single pass of the liquor is usually small 
(loss than 1,0 por cent), consequently the product is recycled 
to achieve a higher concentration of the final product. 
Simulation programme has been run to study the effect of 
recycle ratios of 1.5, 2.5, 3.5 and 4.5 on performance. Feed 
temperature and feed rate has been kept fixed. The effect of 
recycle ratio is shown by tho results in Table 4.12 and 4.13 
for FCui with vertical and horizontal ealendria respectively. 
Figure 4.16 and 4.17 show that product concentration increases 
with increase in recycle ratio. Product concentration 
increases by 12 per cent for a three fold rise in recycle 
ratio. As recycle ratio is increased, more and more of concen- 
trated product is being mixed with fresh feed, and as a result 
final product concentration increases. Total pressure drop 
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shovti a nau with recycle ratio. Prom Table 4. 13 It is 
observed that total pressure drop rise direotly proportional 
with recycle ratio, while for vertical calendria this rise 
is about 65 per cent. This can bo attributed to the fact 
ill units wimi horizontal calendria pressure loss rises duo 
bo increase in end losses, which, varies, as the square of 
b.lu; i- 1 c luid velocity inside the tube. As recycle ratio 
increases from 1.5 to 4.5, liquor velocity inside tube increases 
irom 2.1 bo 4.0 m/s. This rise in velocity greatly increases 
pressure loss at the ends. In evaporators with vertical 
calendria, though pressure loss increases with recycle, some 
pressure is recovered when liquor flows downward and more over 
boiling occurs inside tube for vertical orientation. Boiling 
inside Lube and pressure recovery during downward flow contri- 
bute to the* smaller rise in pressure loss. 

Overall heat transfer coefficient does not change upto 
a recycle ratio of 3.5 beyond which it increases. Inspection 
of t ho values in Tables 4,12 and 4.13 shows a rise in average 
boiling hunt transfer coefficient. Contributing towards 
improvement of overall heat transfer coefficient. 

4.3.4 Effec t of food Liq uor Viscosity ; 

Effect of feed liquor viscosity was studied considering 
bamboo and bagasse black liquors. The viscosity for bagasse 
black liquor at 40 per cent concentration at 110°C is 
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TABLE 4.14 


VISCOSITY OH PEREOfiMAIiCE of eokced 

OliLOni ATI 01 EVAPORATOR ( HOE I ZOITAL 
C ALBUM A) 


Evaporator: L=3m, D=5.05/5.45 cm n = 1 1=3 

Pro oo ou Gondltiona 

System: Bagasse kraft Black Liquor 

P ± » 186 kl/m 2 lijj=0.17 m 2 K/kW T ± = 110° C 
= 40 per cent = 2.5 kg/s/tube 


•D U 

m / s 

5 

cp 

%e 

per §ent 

S o 

percent 

U nb 

¥/m 2 K 

2.5 2.45 

15.57 

10480 

43.45 

45.32 

1545 

4.5 3.53 

20.06 

11927 

46.56 . 

48 . 24 

1700 

6.5 4.48 

29.97 

10270 

49.55 

51.19 

1650 
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8 op. The recycle ratio was varied from 2.5 to 6.5. Results 
12 cp and that for bamboo black liquor at 45 per cent and 
110°C is 8 cp. The recycle ratio was varied from 2.5 to 6.5. 
Results lor bagasse black liquor are summarised in Table 4.14. 


Those results are compared with the values reported in 
L.iii.le; 'i .13, tor bamboo black liquor with a feed concentration 
o I. 4 j pci (a. ii l tor identical values of feed rate and temperature, 
these ro sul b o show that overall heat transfer coefficient is 
jjii.tl l.e i let bagasse black liquor compared to bamboo black 


li'luoi. At; a recycle ratio of 2.5 the mixed feed concentrations 
m e 415.5 and 51.0 per cent for bagasse and bamboo black liquors 
respectively, the corresponding viscosities are 15 and 10 cp 
;it lit) 0. 1'he higher viscosity of bagasse black liquor lowers 
the overall heat transfer coefficient. At a recycle ratio of 4.5, 
the rise in final concentration is 12 and 8 per cent for bamboo 
and bagasse black liquors, respectively. 

Table 4.14 shows a increase in overall heat transfer 
coelfici wilt for bagasse black liquor from 1.54 to 1.70 kW/m 2 k, 
for recycle ratio upto 4.5, and decreases subsequently. 
increases from 10480 to 11920, due to increase in recycle 
ratio. Velocity of liquor velocity inside tube increases 
from 2.45 to 3.50 m/s, when recycle ratio is increased from 
2.5 to 4.5. At the same time viscosity liquor entering the 
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tlio tube increases from 15.5 to 20.0 cp. At a recycle 
ratio ox 6.5, though, velocity increases to 4.5 m/s, viscosity 
of liquor entering the tube rises to 30 cp thereby lowering 
the values of and the overall heat transfer coefficient. 




CHAPTER 5 


SUMMARY 

Ma enema fei cal model developed for steady state 
simulation of evaporation process in this work gives reliable 
prediction ox its performance. The model was applied to 
Badger' s experimental data and the results showed satis— 
factory agreement* The present model predicted boiling 
length and overall heat transfer coefficient in boiling 
zone within 340 }?er cent and the agreement was within ^’50 r €* 
cent for evaporation rate av»4 total heat transferred in boiling 
zono. The model has been applied to simulate performance of 
evaporation in a single vertical tube .operating on bamboo 
kraft black liquor. An increase in feed liquor concentration 
from 35 to 60 per cent decreased the overall heat transfer 
coefficient by 20 per cent. An increase in feed temperature 
( 90-110° 0) raised the overall heat transfer coefficient by 
15 per cent. Two-fold increase in feed rate (5-10 kg/s/tube) 
gave a small increase (10 per cent) in overall heat transfer 
coefficient and pressure drop. Change in concentration of the 
food liquor in a single pass unit was observed to be small (less 
than 1 per cent). The model was also used to study the 
performance of a forced circulation evaporator with recycle. 
Performance for evaporators with horizontal and vertical 
orientation of calendria was studied against variations of 
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APPENDIX A 


Py.Q... PHASE PLOW PRESSURE DROP IN BOILING- ZONE 

Martinelli and Nelson [22] postulated that three 
factors contribute to the pressure drop resulting from the 
two-phase flow of boiling mixture; these factors include the 
o.l .loots of, frictional forces, momentum changes and hydrostatic 

force a, 

(a) trie fciomil Pressure Drop Gradient : 

L o c Idiar t-Mar t inell i [22, 9] have shown that two-phase 
flow behaviour can bo considered under four general categories 
consisting of the following pairs of conditions : 

1. Both liquid and vapor can be turbulent , 1 ( t ) -v ( t ) 

2. Both liquid and vapor can be viscous, l(v) -v(v) 

3. Liquid can be viscous and vapor turbulent l(v)-v(t) 

4. Liquid can be turbulent and the vapor viscous l(t)-v(v) 
Amongst these turbulent flow pattern of both liquid and vapor 
phases will be dominant during forced circulation. For this 
case Lockhart and Martinelli [22] have shown that for isothermal 
two-phase flow the frictional pressure drop can be calculated 
from equation (A-l). 



loo 


*/' 1 TPP “ fictional Pressure drop gradient due 
to two -phase flow 

' AL' g “ ^ ressure dro P gradient only vapor is flowing 

ix P erime *itally determined parameter and is 
correlated by means of the dimensionless 

parameter I. ' 

Alternatively , two-phase frictional pressure drop can also be 

obtained from equation (A-2) 

/ dP \ , , 2 

K AL'TPP “ ^~7Tl'L %tt (A-2) 


( 


/SL J 1 


frictional pressure 'drop gradient if only 


liquid is flowing 
0^ - Parameter related 

i'ho relationship between 0, and 0 

1 g 



^gtt 

y 2-n 
A tt 2 


to dimensionless parameter X. 
is given by equation (A-3) 

(A-3) 


The value of the index-n in equation (A-3) depends upon the 

flow mechanism. 

Dimensionless correlation parameter X.^ obtained by 
dimensional analysis is given by equation (A-4) 


A tt 


( Ik ) 1/(2 - n ) ,H) n/(2 - n) * 

' ri V ( 


_L 

¥. 


(A-4) 


g 


The quality of vapor and liquid mixture is defined by 

equation (A-5) 



101 


x 


W 

__g 

w, +w 

1 g 


(A-5) 


liquation (A-4) can be writt 


en in terms of x as equation (A-6) 


L tt 


ft, V , ) ( — > 


g 


x 


(A-6) 


Index >n* has a value of 0.2-0.25 [21]. A value of 

0.25 in used in the calculation. 

Martinelli-Nelson have modified equation (A-2) to get 

the following equation (A-7). 


/ A f\ 

' AL f IP3? 


(ill) 

k AI j Tj, 


(1-x) 1 * 75 0 2 


ltt 


(A-7) 


liquation (A —7) was obtained as a result of empirically deter- 
mined relationship between superficial liquid pressure drop 
(only liquid flow rate) and the actual liquid pressure drop 

( total flow rate as liquid ) . 


/ A 1’ \ 

' AM/LPF 


(_iU) (i-x ) 1 - 75 

1 ai'i ' ± 1 

o 


(A-8) 


where (-—A-)^ is the frictional pressure drop gradient for 
the flow of 100 per cent liquid. 

The parameter 0 1 ^ has been correlated in terms of 

by Anthony and 'Robert [5]. 


0 itt = Exp[l. 44065 - 0.50445 * In (X tt ) + 0.06212- (ln(X tt ) ' 
-0.00106 • (m(X tt )) 3 - 0.00101- (ln(X tt )) 4 
•K), 00003* (ln(X tt )) 5 + 0.00002* (ln(X tt )) 6 ] 

(A-9) 
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( b ) Accel 0.1/, -i, lion lines 


ssure_ Gradient 


Idle contribution of momentum changes of the fluid ; 
to the accelerative pressure drop dP can be written 


'a 

A*dl> a = Rate of change of momentum 

A«dP. 


as 


' a 


W. tu w u 

&C gf. 


( A-10 ) 


W„ 


W,j,(l-x) 


w g = W T .x 


di ’a onn bu calculated if u ± and u gJ phase velocities are 
known, and u g can be calculated from the knowledge of 

void fraction 'a* . 


W^(l-x) 
1 A(IA a )"p 


u = 


W T x 


g A a ( 


• g 


how substituting for W, , W , u., and u in equation (A-10) we 

JJ g g 

get 


w 2 

dp = 4 - a ( 

a S 0 


(1-x) 2 , X 2 

pTTlta) + TJ 


-) 


s 


/ dP \ 


d /j A=£i4 + _sL) 

u v. p 7 1 -A \ + „ o ' 
5 C Jl 


( A-ll ) 


a g n ” '^(l-a) « f g 

Equation (A-ll) is used to determine accelerative pressure 

drop gradient, 

( c ) Prcutruro Drop Gradient Due to Elevation ; 

Pressure drop gradient due to elevation can be deter- 
mined from the equation; 


/dm 
HFz' 


elevation 


p £ 
•tp g c 


(A-12) 


where j. - weighted average density of two phase flow mixture, 

Jr 
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appendix b 


(a) Pressure drop duo to sudden expansion of vapor-liquid 

mixture. 


long [21] gives equation (B 


/ j- . 


duo to sudden expansion fcr rtea^ater' ^ 




~ -ir 


I) _ D 

*» «j^ 


^Sr« 


2a (1-cr) [ 


a 


JY 

fv 


+ 


%*Ifn 

rr=^r J 


(B-l) 


where a * A l /a 2 

bub scrip Ui 1 and, 2 refer to the positions before and after 
ruotrlo bion, respectively, it is assumed that thermodynamic 
phase librium exist s before and after the restriction and 

that no phase change occurs. 

(b) Ire same d i op duo to sudden expansion and contraction of 


single phase liquid. 


Expansion loss: Pg-P^ 


h U P 

2 g n 


(1-a) 


iA 


Contraction loos: P 0 ~P, = K ~ — 

c X Q c. Qq 

K. q « contraction coefficient 

0.4(1 - ~) for turbulent flow. 


(B-2) 

(B-3) 
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APPENDIX Q 

ANTOINE EQUATION CONSTANTS EQR WATER 


TABLE C-1 


Temperature 
range . 0 0 

— . .. Antoine 

A ~ 

Equation Con 
■ 

stants 

0-30 

8.184254 

-D 

1791.3 

c 

238.10 

30-40 

8.1393986 

1762.262 

236.29 

40-50 

8.0886767 

1739.351 

234.10 

50-60 

8.0464204 

1715.429 

232.14 

60-70 

8.0116295 

1965.167 

230.41 

70-80 

7.9845588 

1678.948 

228 . 97 

80-90 

7.9634288 

1665.924 

227.77 

90-100 

7.948396 

1656.39 

226.86 

100-150 

7.9186968 

1636.909 

224.92 


log 10 P = A ~ 5T+ C 

T = Temper a ture,°C 


P = Pressure, mm of Hg 
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TABLE C-2 


Pressure Range 
mm Hg 


Antoine Eauation 

Constants 

A 

B 

- C 

5-30 

8.184254 

1791.3 

238.1 

30-35 

8.1393886 

1707.262 

236.29 

35-90 

8.0886767 

1739.351 

234.10 

90-150 

8.0464202 

1715.429 

232.14 

150-230 

8.0116295 

1695.167 

230.41 

230-340 

7.9845538 

1678.948 

228.97 

340-520 

7.9634238 

1665.924 

227.77 

520-760 

7.948396 

1656.39 

226.86 

760-1500 

7.9136968 

1636.909 

224.92 


1 " A-log 1Q P 

T sa Temperature, °G 
P = Pressure mm of Hg 
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APPENDIX D 


3MLL SIDE HEAT TRANSFER COE'PFI C IENT 


The shell side heat transfer coefficient for the film 
wiso condensation of steam can be determined from equation D-l 

for lone vortical tube. 


h 


1.13 [ 


«fp ® H q 25 

™ — qjrr-g 

' x arr x a ' 


(D-l) 


" ^f'-gV “s' 

Shell side boat transfer coefficient for filmwise condensation 

on horizontal tube can bo determined from equation D-2. 

■ <m 

a 


K PjS 6 AH 0 25 
h = >‘£ tV-'-'s 1 ] ’ 


(D-2) 


H = latent boat of condensation, BTU/lb 
g - gravitational acceleration ft/hr 2 
L = length of tube, ft 
D q = outside tube diameter, ft 

:= viscosity of condensate film, lb/ft-hr 
T ”• saturated vapor temperature, °P 
» surface temperature, °F 




